Magnetization process of the spin-1/2 XX Z models 
on square and cubic lattices 

Masanori Kohno and Minora Takahashi 

Institute for Solid State Physics, 
University of Tokyo, Roppongi, Minato-ku, Tokyo 106 
(February 7, 2008) 

Abstract 

The magnetization process of the spin-1/2 antiferromagnetic XX Z model 
with Isingdike anisotropy in the ground state is investigated. We show nu- 
merically that the Ising-like XXZ models on square and cubic lattices show 
a first-order phase transition at some critical magnetic field. We estimate 
the value of the critical field and the magnetization jump on the basis of the 
Maxwell construction. The magnetization jump in the Ising-limit is investi- 
gated by means of perturbation theory. Based on our numerical results, we 
briefly discuss the phase diagram of the extended Bose-Hubbard model in the 
hard-core limit. 

PACS numbers: 75.30.Kz, 75.40.Cx, 05.30.Jp, 67.40.Db 
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I. INTRODUCTION 



Neel Q predicted the first-order transition of anisotropic antiferromagnets in the presence 
of a magnetic field in 1936. He pointed out that the spins will abruptly change directions from 
parallel to perpendicular with respect to the c-axis (easy axis of sublattice magnetization) at 
some value of the external magnetic field, when the magnetic field is applied to the direction 
parallel to the c-axis. His prediction was confirmed experimentally 0, and this first-order 
phase transition is now known as the spin-flopping process. Thirty years later than the 
discovery of the spin-flopping process, C.N. Yang and CP. Yang showed by the Bethe ansatz 
that the one-dimensional (Id) spin-1/2 XX Z model with Ising-like anisotropy exhibits a 
second-order transition in the presence of a magnetic field |§. Thus, for one-dimensional 
Ising-like antiferromagnets, quantum fluctuations, which are neglected in the mean-field 
approximation, play an essential role. In this way, quantum fluctuations may drastically 
modify the classical behavior depending on the dimensionality. Hence, we investigate the 
magnetization process of the spin-1/2 Ising-like XX Z (l-XXZ) models in two and three 
dimensions in order to see how quantum fluctuations modify the classical behavior of the 
magnetization process of Ising-like antiferromagnets. 

Also, the spin-1/2 XX Z model can be translated into the hard-core boson model with 
nearest neighbor repulsion JJ[. This model corresponds to a special case of the extended 
Bose-Hubbard model which is considered to be relevant for low-temperature properties of 
liquid helium on a periodic substrate and also for Josephson junction arrays From the 

theoretical point of view, a lot of attention has been paid to the Bose-Hubbard model as 
the simplest model to describe the superfluid-insulator transition |5|-|9|]. We can obtain in- 
formation about the superfluid-insulator transition occurring in the extended Bose-Hubbard 
model through the investigation of the spin-1/2 l-XXZ model. 

This paper is organized as follows: In Sec. ||, the XX Z model is defined. The details of 
the numerical calculations are presented. In Sec. |ITT| , we review the classical Ising-like XX Z 
model. Numerical results on the magnetization curve of the spin-1/2 l-XXZ models in two 
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and three dimensions are shown in Sec. |IV|. In Sec. |V|, the behavior of the magnetization 
curve in the Ising-limit is investigated by means of perturbation theory. In Sec. [VT] , we 
briefly discuss the superfluid-insulator transition in the extended Bose-Hubbard model in 
the hard-core limit based on the numerical results of the spin-1/2 1-XXZ model. Section 
Vll| is devoted to summary. 



II. MODEL AND METHOD 

In the present paper, we consider the XX Z model defined by the following Hamiltonian: 
'Hxxz = J + SfSj + (2-1) 

(id) 

where denote the x(y, z) components of the spin operator at site i. Here denotes 

nearest neighbors. The anisotropic coupling constant is denoted by A. For A = 1, the 
isotropic Heisenberg model is recovered. We investigate the spin-1/2 XX Z models on square 
and cubic lattices in the ground state in the canonical ensemble. Namely, we measure the 
energy E within the subspace of fixed magnetization M (=J2i S*). The magnetic field in a 
finite-size cluster is defined as H(M) = {E(M 1 ) - E(M 2 ))/{M l - M 2 ), where M is defined 
by (Mi + M 2 )/2. In the thermodynamic limit, this definition of the magnetic field reduces 
to the normal one: H = dE/dM. The maximum magnetization and the saturation field 
are denoted by M max (=iV s /2) and iJ max (=J(A + l)d), respectively. Here N s and d are the 
system size and the spatial dimensionality. 

We use the Lanczos algorithm (exact diagonalization) for clusters up to 32 sites and 
the cluster algorithm EDI (quantum Monte Carlo) for larger clusters up to 100 sites. For 
the cluster algorithm, the measurements have been performed at the inverse temperature 
(3 J = 16. The width of the Trotter slice is chosen as At J = 0.04 for two dimensions (2d) 
and At J = 0.053 for three dimensions (3d). The simulation has been performed in the 
canonical ensemble. In the small A regime, the number of points near M = obtained by 
exact diagonalization is too small to use the Maxwell construction. On the other hand, in 
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the large A regime, statistical errors in the quantum Monte Carlo calculation become large, 
because the spin configurations are almost frozen. Hence, we have investigated the small A 
regime by quantum Monte Carlo and the large A regime by exact diagonalization. 

In order to see finite-size effects, we show the size-dependence of the energy gap A g (= 
EiM = 1) — E(M = 0)), the critical field H c and the magnetization jump M s in Fig.[TJ. As 
shown in this figure, the size-dependence is very small. As a check of numerical accuracy, we 
compare the energy gap A g obtained by this method with those obtained by other methods 
as shown in Fig.|2|. Our result is quite consistent with those of third-order spin-wave theory 



TT| and the series expansion around the Ising-limit ||12|| . Hence, we consider that the inverse 
temperature (3 and the width of the Trotter slice At are sufficient. [See also Fig.f|.] For 
three dimensions, we report the 64-site results. 

III. REVIEW OF THE CLASSICAL SPIN CASE 

Before investigating the spin-1/2 XX Z models, we briefly review the magnetization 
process of the classical l-XXZ model in the ground state [0,0]. The ground state energy 
of the classical l-XXZ model may be written in the following form: 

Eq 

' /X : ' S \sin(# + 0) sin(0-0) 
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+Acos(# + 0)cos(#-0)), (3.1) 

where S^fej^ represent the x(y, z) components of the spin operators at a site in the A(B) 
sublattices. The length of the spin and the coordination number are denoted by S and z(= 
2d), respectively. The angles 9 and are defined as in Fig||(a) (0 < 6 < n/2, < <p < vr/2). 
The Zeeman term E% is written in the following form: 

E z = ~^(S Z A + S* B ) = ^^(cos(6 + 0) - cos(9 - 0)). (3.2) 

By minimizing the total energy (E tot = Eq-xxz + Ez) with respect to 9 and 0, one finds 
the following stable states (Fig|J(b)): 



(i) e = 0, = 0, E tot = -JX for H < J^fW^l 

(ii) 9 = tt/2, = arcsin -J^- y E tot = —J — -Jg^ for J(\ + 1)>H> J^fW^l 

(iii) = tt/2, = tt/2, E tot = JX-2H for if > J(A + 1), 

where J and if are defined as J = JN s zS 2 /2 and if = HN s S/2. The magnetization curve 
of the classical 1-XXZ model is shown in Fig.^|(c). The transition from the state (i) to the 
state (ii) is known as the spin-flopping process P|,|T3||. The critical field H c is defined as the 
magnetic field above which the ground state has non-zero magnetization. The A dependence 
of the critical field H c and that of the magnetization jump M s are obtained as 

H c /H max = M s /M max = ^/(A-1)/(A + 1), (3.3) 

where if max = JSz(X + 1) and M max = N S S. 



IV. MAGNETIZATION CURVE OF THE SPIN-1/2 XXZ MODEL 

In this section, we present the numerical results on the magnetization curve of the spin- 
1/2 1-XXZ models on square and cubic lattices. As an example, we show the magnetization 
curve of the spin-1/2 XXZ model for A = 2 on a square lattice in Fig. The critical field 
H c and the magnetization jump M s are determined on the basis of the Maxwell construction 
as follows. [See also Fig.||(a).] We fit the energy as a function of magnetization by a 
polynomial. The tangent from the point at M = to the fitting curve gives a lower energy 
than the numerical data in the region of < M < M s . Here M s is the magnetization at 
the point of contact between the fitting curve and the tangent. Hence, we can identify the 
region of phase separation as < M < M s on the basis of the Maxwell construction |I4 



The magnetic field of the phase-separated state (H c ) is given as the slope of the tangent. In 
practice, we have determined the phase-separation boundary as the point where the following 
condition is satisfied: dE/dM = (E(M) — E(M = 0))/M as in Fig.[5]. The size-dependence 
of the critical field H c and that of the magnetization jump M s determined by the Maxwell 
construction are very small as discussed in Sec.[TT] (Fig .[I]). 
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We show the A dependence of the critical field H c in Fig.||. The critical field H c is sup- 
pressed by quantum fluctuations. In order to see how large the critical field H c is suppressed 
by quantum fluctuations, we have tried to fit the numerical data as H c /H max = (jtj) 01 , 
analogously with the classical result a = 0.5 (eq.(pT3|)). We estimate a to be a = 0.64±0.01 
for 2d and a = 0.57 ± 0.01 for 3d. Note that the A dependence of the critical field H c in 
two and three dimensions is quite different from the one-dimensional case, where the gap 



(—H c = dE/dM\ M / Mm!i ^ + o) opens exponentially: H c oc exp[— 7r 2 /2y2(A — 1)] |L5 |. 

Here, we mention the relation between the critical field H c and the energy gap A g . It is 
expected that the energy gap A g is larger than the critical field H c , if a first-order transition 
occurs in the presence of a magnetic field. The reason is as follows. The ground state of 
M = 1 is considered to be the one-magnon state, which may be described by spin-wave 
theory. Hence, the gap A g corresponds to the excitation energy of one magnon from the 
ground state of M = 0. On the other hand, phase separation occurs, because magnons gain 
energy by interacting attractively with each other. The critical field H c would be determined 
by the effective attractive interactions between the macroscopic number of magnons. As a 
result, if phase separation occurs, the gap A g is expected to be larger than the critical field 
H c , i.e. 

A g > H c = dE/dM\ M/Mm ^ +0 , (4.1) 

where M is assumed to be a macroscopic number when the limit M/M max — > +0 is taken. 
We compare the gap A g and the critical field H c of the spin-1/2 I-XXZ model on a square 
lattice in Fig.|7|. The gap A g is always larger than the critical field H c as expected. It is 
interesting to contrast this behavior with the one-dimensional result. For one dimension, the 
transition is of second order JJ, and the following relation is satisfied: dE /dM\M/M max -^+o = 
E(M = 1) — EiM = 0). This is considered to be due to effective repulsive interactions. 

The A dependence of the magnetization jump M s is shown in Fig.|8|. We estimate the 
critical value of A, where M s vanishes, as A c = 1.00 ± 0.02 by extrapolating the data in 
Fig||. This confirms that the spin-1/2 l-XXZ models on square and cubic lattices show a 
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first-order transition at some critical field for any value of the anisotropic coupling constant 
larger than one (A > 1). The A dependence of the magnetization jump M s is remarkably 
different from the classical result, especially in the large A regime. 



V. ISING-LIMIT 

In this section, we discuss the magnetization process of the spin-1/2 XX Z model in the 
Ising-limit. In Fig|| the value of the magnetization jump M s in the Ising-limit (A — > oo) 
does not coincide with that of the Ising model (A = oo)(M s (A = oo) = M max ): 

M S (A -> oo) ^ M S (A = oo). (5.1) 

This can be explained by means of perturbation theory as follows. We rewrite the spin-1/2 
XX Z model as 

Hxxz = ./»<- s i + ' J ^ S '' S 'J + S ? S V- ( 5 - 2 ) 

where J and e are defined as J = JA and e = 1/A. We consider the XF-term as the 
perturbation. At M = 0, the unperturbed ground states are the two-degenerate Neel states. 
The leading perturbation energy is of order e 2 . On the other hand, in the limit of M — > M max , 
the leading perturbation energy is of order e and proportional to M — M max . Hence it 
is expected that phase separation occurs for the magnetization smaller than some value 
M s (< M max ) in the Ising-limit. We numerically estimate the value of M s in the Ising-limit 
with first-order perturbation theory in the following way. The first-order perturbation energy 
Ei is obtained as 

„ eJJ2 a> p(P\ E(i,j){S?S- + S?S]~)\a) 

&i - -7; / i \ ' y b - 6 

2 22a\ a \ a ) 

where \a) and \(3) denote unperturbed ground states of the Ising model in the subspace 
of fixed magnetization. We generate |a)'s randomly and measure Ei using Monte Carlo 
technique. The value of M s is determined based on the Maxwell construction. Figure |9| 
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shows the first-order perturbation energy E\ and the value of M s in the Ising-limit on hyper- 
cubic lattices in dimensions up to six. We extrapolate the data in Fig.|9|(b) and estimate the 
inverse dimensionality, where M s coincides with M max , as 1/d = 0.01 ± 0.02. This confirms 
that the value of M s in the Ising-limit (A — > do) does not coincide with that of the Ising 
model (A = do) in finite spatial dimensions. 

VI. RELATION TO THE BOSE-HUBBARD MODEL 

In this section, we briefly discuss the phase diagram of the hard-core boson model with 
nearest-neighbor repulsion based on the numerical results of the spin-1/2 XX Z models. The 
hard-core boson model with nearest-neighbor repulsion can be obtained from the following 
extended Bose-Hubbard model by taking the on-site repulsion U infinity: 

Hbh = t J2(blbj + b]k) + Uj2n i (n i - 1) + V ^(m - 1/2) (n,- - 1/2), (6.1) 

where b\ (&») creates (annihilates) a boson on site % and rij = b\bi. Here denotes nearest 
neighbors. The spin-1/2 XX Z model can be translated into the hard-core boson model 
with nearest-neighbor repulsion (t <-» J/2, V <-> JA/4) [|J. Figure |8] corresponds to the 
phase diagram of this model by relating the filling p and the interaction strength V/t to 
M/M max and A according to p — (1 — M/M max )/2 and V/t = A/2. The numerical results 
in the previous section are translated as follows: The superfluid-insulator transition occurs 
in the region of V > t/2. This transition is a first-order transition, which is consistent with 
recent investigation of the Bose-Hubbard model |5) . Phase separation does not occur for the 
density p smaller than p c , for finite V/t, where p c = (1 — M S (A — > oo)/M max )/2 > 0. This 
p c approaches zero as the spatial dimensionality d goes to infinity (Fig.^(b)). 

VII. SUMMARY 

In summary, numerical results on the magnetization process of the spin-1/2 Ising-like 
XX Z models have been reported. The spin-1/2 XX Z models on square and cubic lattices 
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show a first-order phase transition at some critical magnetic field for the anisotropic coupling 
constant larger than one (A > 1). The critical field H c and the magnetization jump M s are 
estimated on the basis of the Maxwell construction. The critical field H c is suppressed by 
quantum fluctuations (Fig.|6|). We have demonstrated that the energy gap A g is larger than 
the critical field H c (Fig.0). The anisotropy A dependence of the magnetization jump M s is 
remarkably different from the classical result (Fig||). It is strongly suggested that the value 
of M s in the Ising-limit (A — > oo) does not coincide with that of the Ising model (A = oo) in 
finite spatial dimensions due to quantum effects (Fig.|). 
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FIGURES 

FIG. 1. Size dependence of the critical field H c [(a)], the magnetization jump M s [(b)] and the 
energy gap [(c)] A g on the two-dimensional spin-1/2 XX Z model. Dotted lines are guides to the 
eye. 

FIG. 2. Energy gap A g for the two-dimensional spin-1/2 XX Z model. Open circles and open 
squares denote the results obtained by third-order spin-wave theory cited from [jllj and the series 
expansion cited from (T^j, respectively. Solid and dotted lines correspond to the classical and 
one-dimensional cases, respectively. 

FIG. 3. (a) Definition of 9 and <ft. (b) Schematic picture of stable spin-configurations, (i) Neel 
state, (ii) spin flopping state, (iii) fully polarized state, (c) Magnetization curve of the classical 
XX Z model with Ising-like anisotropy. 

FIG. 4. Magnetization curve of the two-dimensional spin-1/2 XX Z model with Ising-like 
anisotropy. Open and solid symbols denote the data in 26-site and 36-site clusters, respectively. 
Solid (dotted) lines correspond to the phase-separated states determined on the basis of the Maxwell 
construction by using the data in 36-site (26-site) clusters. We choose M\ — M2 = 1 or 2 in the 
definition of the magnetic field H. The definition of the magnetic field H is in the text. The 36-site 
data for < M/M max < 0.6 are obtained by quantum Monte Carlo. Other data are obtained by 
exact diagonalization. 

FIG. 5. dE/dM (solid symbols) and (E(M)-E(M = 0))/M (open symbols). The intersection 
point of two curves corresponds to the phase-separation boundary. The data are obtained by exact 
diagonalization in a 26-site cluster (A = 2). We set J = 1. 

FIG. 6. Anisotropy A dependence of the critical field H c of the spin-1/2 XX Z model, (a) Linear 
plot and (b) Log-log plot. Solid and dotted lines correspond to the classical and one-dimensional 
cases, respectively. Dashed and dash-dotted lines correspond to the cases a = 0.64 and a = 0.57. 
The definition of a is in the text. 
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FIG. 7. Anisotropy A dependence of the energy gap A g (solid symbols) and the critical field H c 
(open symbols) of the spin-1/2 XX Z model on a square lattice. Solid and dotted lines correspond 
to the classical and one-dimensional cases, respectively. 

FIG. 8. Anisotropy A dependence of the magnetization jump M s of the spin-1/2 XX Z model. 
The solid line corresponds to the classical case. Bold lines are guides to the eye. This figure 
corresponds to the phase diagram of the hard-core boson model with nearest-neighbor repul- 
sion by relating the filling p and the interaction strength V/t to M/M max and A according to 
p = (1 — M/M max )/2 and V/t = A/2. The definition of the hard-core boson model is in the text. 

FIG. 9. (a) First-order perturbation energy in the Ising-limit on hyper-cubic lattices. We set 
J = 1 as the energy unit. The system sizes are 400 sites for 2d, 1000 sites for 3d, 1296 sites for 
4d, 1024 sites for 5d and 4096 sites for 6d. Dashed lines denote the tangents from M = 0. The 
magnetization at the point of contact corresponds to M s . (b) Dimensionality d dependence of M s 
in the Ising-limit on hyper-cubic lattices. The dashed line represents M s /M max = 1/d. 
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